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(57) Abstract: Methods and apparatus compress data, comprising an In-phase (I) component and a Quadrature (Q) component. The 
compressed data (212) may be saved into a memory or may be transmitted to a remote location for subsequent processing or storage. 
Statistical characteristics (505, 507) of the data are utilized to convert the data into a form that requires a reduced number of bits in 
accordance with the statistical characteristics. The data may be further compressed by transforming the data, as with a discrete cosine 
transform (203), and by modifying the transformed data in accordance with a quantization conversion table (205) that is selected 
^ using a data type associated with the data. Additionally, a degree of redundancy may be removed fix)m the processed data with an 
^ encoder (207). Subsequent processing (1200) of the compressed data may decompress the compressed data in order to approximate 
1^ the original data by reversing the process for compressing the data with corresponding inverse operations. 
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MEASUREMENT AND SIGNATURE INTELLIGENCE 
ANALYSIS Am REDUCTION TECHNIQUE 

This application claims priority to U.S. Provisional Application Serial. No. 60/392,316 
entitled 'Measurement and Signature Intelligence Analysis and Reduction Technique', and filed 
June 28, 2002, and U.S. Application No. 10/269,818, entitled 'Measurement and Signature 
Intelligence Analysis and Reduction Technique \ and filed October 11, 2002. 

FIELD OF THE INVENTION 

The present invention relates to compressing and decompressing data such as synthetic 
aperture radar data. 

BACKGROUND OF THE INVENTION 

Compression of Synthetic Aperture Radar (SAR) data may require that both magnitude 
and phase information be preserved. Figure I shows data processing of syntiietic aperture radar 
data according to prior art. Synthetic aperture radar data 102 are typically collected in analog 
fonnat by an antenna 101 and is converted to digital format through an Analog-to-Digital (A/D) 
converter 103. The raw, unprocessed data af e refenred to as Video Phase History (VPH) data 
104, and comprise two components: In-phase CO and Quadrature (Q). Video phase history data 

104 having multiple components, such as I and Q, are typically referred as complex SAR data. 
Complex SAR data are essential for the generation of complex SAR applications products such 
as interferograms, polarimetry, and coherent change detection, in which a plurality of such 
images must be processed and compared. 

Video phase history data 104 are then passed through a Phase History Processor (PHP) 

105 where data 104 are focused in both range (corresponding to a range focusing apparatus 107) 
and azimutii (corresponding to an azimuth focusing apparatus 109). The output of phase history 
processor 105 is referred to as Single Look Complex (SLC) data 110. A detection fimction 111 
processes SLC data 1 10 to form a detected unage 1 12. 

Existing complex SAR sensors collect increasingly large amounts of data. Processing 
the complex data information and generating resultant imagery products may utilize four to eight 
times the memory storage and bandwidth that is required for the detected data (I&Q). In fact. 
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some studies suggest exponential growth in associated data throughput over the next decade. 
However, sensors are typically associated with on-board processors that have limited processing 
and storage capabilities. Moreover, collected data are often transmitted to ground stations over a 
radio channel having a limited frequency bandwidth. Consequently, collected data may require 
5 compression in order to store or transmit collected data within resource capabilities of data 
collecting apparatus. Also, a SAR compression algorithm should be robust enough to compress 
both VPH data 104 and SLC SAR data 110, should produce visually near-lossless magnitude 
image, and shoxild cause minimal degradation in resultant products 1 12. 

Several compression algorithms have been proposed to compress SAR data. However, 
10 while such compression algorithms generally work quite well for magnitude imagery, the 
compression algorithms may not efficiently compress phase information. Moreover, the phase 
component may be more important in carrying information about a SAR signal than the 
magnitude component With SAR data 102, compression algorithms typically do not achieve 
compression ratios of more than ten to one without significant degradation of the phase 
15 information. Because many of the compression algorithms are typically designed for 
Electro/Optical (EO) imagery, the compression algorithms rely on high local data correlation to 
achieve good compression results and typically discard phase data prior to compression. Table 1 
lists several compression algorithms discussed in the literature and provides a brief description 
of each. 



Table 1 : Popular Alternative SAR Data Compression Algorithms 


Compression Algorithm 


Description 


Block Adaptive Quantization (BAQ) 


Choice of onboard data compression 
methods due to simplicity in coding and 
decoding hardware. Low compression 
ratios achieved (< 4: 1). 


Vector Quantization (VQ) 


Codebook created assigning a number for a 
sequence of pixels. Awkward 
implementation since considerable 
complexity required in codebook 
formulation. 


Block Adaptive Vector Quantization 
(BAVQ) 


Consists of first compressing data with 
BAQ and then following up with VQ. 
Similar to BAQ. 


Karhunen-Loeve Transform (KLT) 


Statistically optimal transform for 
providing imcorrelated coefficients; 
however, computational cost is large. 
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Table 1 : Poptilar Alternative SAR Data Conq?ression Algorithms 


Compression Algorithm 


Description 


Fast Fourier Transform BAQ (FFT- 
BAQ) 


2-D Fast Fourier Transform (FFT) 
performed on raw SAR data. Before raw 
data is transformed, dynamic range for 
each block is decreased using a BAQ. 


Uniform Sampled Quantization (USQ) 


Emphasizes phase accuracy of selected 
points. 


Flexible BAQ (FBAQ) 


Based on minimizing mean square error 
between original and reconstructed data. 


Trellis-Coded Quantization (TCQ) 


Unique quantizer optimization design. 
Techniques provide superior signal to noise 
ratio (SNR) performance to BAQ and VQ 
for SAR. 


Block Adaptive Scalar Quantization 
(BSAQ) 


BSAQ's adaptive technique provides some 
perforaaance improvement. 



Existing optical algorithms are inadequate for compressing complex multi-dimensional 
data, such as SAR data compression. For example with optical imagery, because of a human 
eyesight*s natural high frequency roU-ofif, the high frequencies play a less important role than 
low frequencies. Also, optical imagery has high local correlation and the magnitude component 
is typically more important than the phase component However, such characteristics may not be 
applicable to complex multi-dimensional data. Consequently," a method and apparatus that 
provides a large degree of compression without a significant degradation of tihe processed signal 
are beneficial in advancing the art in storing and transmitting complex multi-dimensional data. 
Furtheraaore, the quality of the processed complex multi-dimensional data is not typically 
visually assessable. Thus, a means for evaluating the effects of compression on the resultmg 
processed signal is beneficial to adjusting and to evaluating the compression process. 

BRIEF SUMMARY OF THE INVENTION 

The present invention provides methods and apparatus for compressing data comprising 
an In-phase (£) component and a Quadrature (Q) component. The compressed data may be 
saved into a memory or may be transmitted to a remote location for subsequent processing or 
storage. Statistical characteristics of the data are utihzed to convert the data into a form that 
requires a reduced number of bits in accordance with its statistical characteristics. The data may 
be further compressed by transforming the data, as with a discrete cosine transform, and by 
modifying the transformed data in accordance with a quantization conversion table that is 



wo 2004/004309 



PCTAJS2003/020295 



4 

selected using a data type associated with the data. Additionally, a degree of redundancy may be 
removed from the processed data with an encoder. Subsequent processing of the compressed 
data may decompress the compressed data in order to approximate the original data by reversing 
the process for compressing the data with corresponding inverse operations. 

5 In a first embodiment of the mvention, data are compressed with an apparatus 

comprising a preprocessor, a transform module, a quantizer, an encoder, and a post-processor. 
The preprocessor separates the data into an I component and a Q component and bins each 
component according to statistical characteristics of the data. The transform module transforms 
the processed data into a discrete cosine transform that is quantized by the quantizer using a 

10 selected quantization conversion table. The encoder partially removes redundancy from the 
output of the quantizer using HufBnan coding. The resulting data can be formatted by a post- 
processor for storage or transmittal. With a second embodiment, the preprocessor converts the I 
and Q components into amplitude and phase components and forms converted I and Q 
components. 

15 Variations of the embodiment may use a subset of the apparatus modules of the first or 

the second embodiment In a variation of the embodiment, the apparatus comprises a 
preprocessor, a transform module, and a quantizer. 

BRIEF DESCRIPTION OF THE DRAWINGS 

A more complete understanding of the present invention and the advantages thereof may 
20 be acquired by referring to the following description in consideration of the accompanying 
drawings, in which like reference numbers indicate like features and wherein: 

Figure 1 shows data processing of synthetic aperture radar data according to prior art; 

Figure 2 shows an apparatus for compressing data in accordance with an embodiment of 
the invention; 

25 Figure 3 shows a preprocessor apparatus for preprocessing a complex image in 

accordance with an embodiment of the invention; 

Figure 4A shows a process for binning data associated witb a complex image in 
accordance with an embodiment of the invention; 
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Figure 4B shows a process for truncating magnitude and phase components of a complex 
image in accordance with an embpdiment of tiie invention; 

Figure 5 shows probability density functions that are associated with In-phase (I) and 
Quadrature (Q) components of exemplary synthetic aperture radar (S AR) data; 

5 Figure 6 shows Root Mean Square Error (RMSE) values that are associated with 

magnitude and phase data for processed signal data as shown in Figure 2 in accordance with an 
embodiment of the invention; 

Figure 7 shows a partitioning of complex image data in order to obtain Discrete Cosine 
Transform (DCT) in accordance with an embodiment of the invmtion; 

10 Figure 8 shows an apparatus for quantizing Discrete Cosine Transform (DCT) data in 

accordance with an embodiment of the invention; 

Figure 9 shows a representative histogram for a low order Discrete Cosine Transform 
(D^T) coefficient in accordance with an embodiment of the invention; 

Figure 10 shows a representative histogram for a high order Discrete Cosine Transform 
1 5 (DCT) coefficient in accordance with an embodiment of the invention; 

Figure 1 1 shows a heuristic process for determining a quantization matrix according to 
an embodiment of the invention; and 

Figure 12 shows an apparatus for decompressing data in accordance with an embodiment 
of the invention. 

20 DETAILED DESCRIPTION OF THE INVENTION 

In the following description of tiie various embodiments, reference is made to the 
accompanying drawings which form a part hereof, and m which is shown by way of illustration 
various embodiments in which the invention may be practiced It is to be understood that other 
embodiments may be utilized and structural and functional modifications may be made without 
25 departing from the scope of the present invention. 

Figure 2 shows an apparatus 200 for compressing Synthetic Aperture Radar (SAR) data 
202 in accordance with an embodunent of the invention. Synthetic Aperture Radar (SAR) data 
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202 can be compressed by apparatus 200 from Video Phase History (VPH) data format 104 or 
from a processed version typically referred as a single look complex SLC fonnat 110. There are 
advantages and disadvantages associated with each format VPH data 104 is available almost 
immediately, but is highly uncorrelateA Single look complex SLC data 110 exhibits some local 
correlation. SLC data 110 may yield slightly better compression results than with VPH data 104, 
but SLC 1 10 data are only available after processing has occurred. 

Other embodiments of the invention may support other applications of complex 
multidimensional data, including weather data, oil and gas exploration data, encrypted/decrypted 
data, medical archival of MRI/CTI and tiiree dimensional sonograms, digital video signals, and 
modem appUcations. 

Referring again to Figure 2, apparatus 200 comprises a preprocessor 201, a transform 
module 203, a quantizer 205, an encoder 207, and a post-processor 209 in order to provide 
compressed data 212. SAR data 202 may comprise SAR pixel data that may be provided in the 
form of two floating-point numbers representing In-phase (I) and Quadrature (Q) components. 
(SAR data 202 may be considered as being "received" even though the data may not be received 
from a radio receiver but obtained from a memory that stores the data.) Preprocessor 201 may 
convert the I and Q components to Magnitude (M) and Phase (cp) components in accordance 
with a second embodiment as will be discussed in the context of Figure 4B. Additionally, 
preprocessor 201 may convert the I and Q components into magnitude and phase components to 
facilitate viewmg SAR data 202. The magnitude and tiie phase components may be obtained 
from the in-phase and quadrature components by using Equations 1 and 2. 

m = (iHqY^ (EQ. 1) 
9-tan"^(Q/I) (EQ.2) 

Moreover, I and Q components may be obtained fit)m the magnitude and the phase 
components by using Equations 3 and 4. 

I = M cos 9 (EQ. 3) 
Q = Msm<p(EQ.4) 

Additionally, the power of a SAR signal may provide good visual results when printing 
intensity (magnitude-only) imagery. The power of a SAR signal may be obtained from Equation 
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P = 201ogioM^(EQ. 5) 

The conversion between (I, Q) and (M, <p) as expressed in EQs. 1-4 allows SAR data 202 
to be studied in both data formats before and after compression. When SAR data 202 are 
represented as magnitude and phase componmts, additional bits may be allocated to the phase 

5 component versus the magnitude component to achieve the least degradation of the phase 
product, depending on characteristics of SAR data 202. In an embodiment, more bits (e.g. as six 
bits) of the phase component and fewer bits (e.g. two bits) of the magnitude component are used 
to generate compressed I and Q components. Conversely, when SAR data 202 are represented 
by in-phase and quadrature components, apparatus 200 can process the in-phase component 

10 separately from the quadrature component for a single complex image. 

Preprocessor 201 also determines a data type (as discussed in the context of Figure 3) 
and informs quantizer 205 through an adaptive control loop 251. 

Figure 3 shows preprocessor apparatus 201 (as shown in Figure 2) for preprocessing 
complex image 202 in accordance with an embodiment of the invention. Preprocessor apparatus 

15 201 reduces the number of bits that are needed to represent complex data (I,Q) within a 
specified degradation (corresponding to an error metric). Typically, VPH data 104 or SLC 110 
data are represented by (I,Q) data pairs 202, in which each pair uses 64, 32, or 16 bits, and where 
I and Q are separately represented in 32, 16, . or 8-bit formats, respectively. Data 202 may be 
formatted in which an ordering of the most significant to the least significant bytes may be 

20 reversed with respect to the assumed order that preprocessor 201 processes data 202. In such a 
case, preprocessor 201 may perform "byte swapping" to reorder data 202 in accordance with the 
assumed ordering of the constituent bytes. 

An adaptive source data calculations module 301 separately processes the I and Q 
components of (I,Q) data pairs 202 in order to determine corresponding statistical 

25 characteristics. (An example of statistical characteristics is shown in Figure 5, in which the I 
component has approximately the same statistical charact^stics as the Q component.) In the 
embodiment, a general-purpose computer (e.g. an associated microprocessor) measures the 
number of occurrences of the I component or the Q component as a Amotion of the value of the I 
component or the Q component Additionally, adaptive source data calculations module 301 

30 performs header analysis by reading information provided at the beginning of a data file 
comprising data 202 in order to determine the format of the data being analyzed, e.g. the number 
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of bits that are associated with (I,Q) data 202. Module 301 also performs data analysis that 
provides statistical characteristics of data 202 as may be characterized by probability density 
functions of the I component and the Q component (as exemplified by Figure 5). Module 301 
determines a bin assignment that may vary with the a value of the I or Q component. In the 

5 embodiment, a size of a bin is inversely related to a value of the probability density function at a 
midpoint of the bin. A calculations module 303 uses the statistical characteristics of data 202 to 
assign the I and Q components into bins. A module 305 uses the bin identity to form tiie V and 
Q' components (converted I component and converted Q component, corresponding to data 204 
in Figure 2), having 8-bit integer values between 0 and 255 by efficiently allocating bins, in 

10 which most of the bins are assigned to a range containing the most data points. For example, 
data (corresponding to either I or Q) may range firom -10,000 to +10,000 units, in which over 
99.9% of the data are contained with a range of -2,000 to +2,000 units. In such a case, most of 
the bins would be allocated between the smaller range (i.e. -2000 to +2,000 units) rather than 
the larger range (i.e. -10,000 to +10,000 units). 

15 a variation of the embodiment, Single Look Complex (SLC) data 1 10 are transformed 

using a Fast Fourier Transform (FFT) prior to binning data 202 by modules 303 and 305, 
wherein a transformation of SLC data 110 has statistical characteristics that are similar tso VPH 
data 104. Qn the embodiment, modules 303 and 305 bin data 202 by first processing the I 
component and subsequently processing the Q component.) However, other embodiments of the 

20 invention may utilize other transform types in order to modify statistical characteristics of the 
data. After quantization by modules 303 and 305, the transformed SLC data are inversely 
transformed using an Inverse Fast Fourier Transform (DFFT). 

Figure 4A shows a process for binning data associated with complex image data 202, as 
performed by module 303 in accordance with an embodiment of the invention. Complex image 

25 data 202 are separated into I and Q conq)onents by a module 403. If the most to the least 
significant bytes need to be reordered, modules 405 and 413 swap bytes for the I component and 
Q component, respectively. Modules 407 and 415 detennine the probability density fianctions 
for the I component and the Q component, respectively over data files (comprising static images 
of a data gathering session). As discussed in the context of Figure 5, the probability density 

30 functions of the I component and the Q component may be essentially the same so that 
embodiments of the invention may utilize one module by separately processmg the I and Q 
components. Modules 409 and 417 bin the I component and the Q component, respectively. The 
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greater the probability density function p(xi), where Xi is the center value of the i'^ bin, the 
smaller the range of the i"" bin in order to provide better resolution for data within the i^ bin. 

Figure 4B shows a process for truncating magnitude and phase components of a complex 
image in accordance with a second embodiment of the invention. In a second embodiment of the 

5 invention, module 305 of preprocessor apparatus 201 may utilize a different number of bits that 
are associated with the phase component (q)) than is associated with the magnitude component 
(M). In the embodiment, fewer bits from the magnitude component (a truncation of M) and more 
bits from the phase component (a truncation of 9), as determined from Equations 1 and 2 by 
converting I and Q into M and 9, are used to generate compressed components V and Q', as 

10 determined from Equations 3 and 4 by converting the truncations of M and 9 into V and Q'. 
Allocating more bits from the phase component helps preserve phase information, as may be the 
case with Video Phase History (VPH) data 104. As shown in Figure 4B, complex image data 
202 is separated into I and Q components by module 453. The I and Q components are 
converted into magnitude and phase components by module 455. Module 457 truncates the 

15 magnitude and phase components in order to retain a desired number of bits from each of the 
components. Module 459 converts the truncated portions of the magnitude and phase 
components to form compressed components P and (corresponding to data 461). 

Apparatus 200 may use the same statistical modeling for the In-phase (I) and Quadrature 
(Q) components if both components have approximately the same statistical characteristics. 

20 Figure 5 shows probability density functions that are associated with in-phase and quadrature 
components of exemplary synthetic aperture radar data. A number of pixels 501 is shown in 
relation to a corresponding pixel values 503 for a typical SAR image. A Probability Density 
Function (PDF) 507 for the in-phase component and a probability density fimction 505 for the 
quadrature component are approximately the same. Figure 5 suggests that apparatus 200 may 

25 process both the in-phase component and the quadrature components in the same way without 
incurring a large error. If probability density fimction 507 is essentially the same as probability 
density function 505, then one may obtain a probability density of one of the components (either 
PDF 507 or PDF 505) and approximate the probability density fimction of the other component 
by the obtained probability density ftinction. However, other embodiments of the invention may 

30 use different statistical relationships for the in-phase component and the quadrature component 
if the statistics characteristics differ appreciably. 
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Prq>rocessor 201 accommodates different sensor types regarding a data format and a 
number of bits per pixel. (A pbcel corresponds to a point in the corresponding image being 
scanned by a radar system.) SAR data 202 are typically 64 bits (with 32 bits for the I component 
and 32 bits for tiie Q component for each pixel) or 32 bits (with 16 bits for the I component and 
16 bits for the Q component for each pixel). Preprocessor 201 determines the range of pixel 
values and the best bin assignment Values of the I and Q components are converted to 8-bit 
formats with more bits being allocated fiom the associated phase component than the magnitude 
component before reducing the I and Q components to 8-bit formats. (As previously discussed, 
two bits from the magnitude component and six bits from tibie phase component, as determined 
from Equations 1 and 2 by converting I and Q into M and cp, are used to generate compressed 
components T and Q', as determined from Equations 3 and 4 by converting the truncations of M 
and cp into r andQ'.) 

Figure 6 shows Root Mean Square Error (RMSE) values that are associated with 
magnitude and phase data for processed data (e.g. processed SAR data 204) as shown in Figure 
2 in accordance with an embodiment of the invention. (The root mean square error is a measure 
of the quantization error by relating the compressed data with the original data.) Values 601 are 
related to an assigned number of bits per pixel 603 for phase data 605, magnitude data 607 (with 
a linear-log representation), and magnitude data 609 (no linear-log representation). Similarly, 
calculations may be performed for (I, Q) data. Root mean square error and Peak Signal to Noise 
Ratio (PSNR) figures of merit may be initially used as a basis for designing preprocessor 201 
and for the evaluating the compressed imagery. 

Processed SAR data 204 (comprising a converted I component and a converted Q 
component) are further processed through transform module 203 using a Discrete Cosine 
Transform (DCT) in order to obtain the frequency representation of the in-phase and the 
quadrature data as transformed data 206 (comprising a transformed I component and a 
transformed Q component). As will be discussed in the context of Figure 7, the converted I 
component and the converted Q component of SAR data 204 are separately partitioned into 
smaller blocks. (Each block is essentially independent of other blocks so that each block may be 
processed individually in order to process an entire image.) The discrete cosine transform is 
well known in the art, and is given by Equation 6. 
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Ni - 1 - 1 



— + 1) 
2-Ni 



cos 



(2 j ^ 1) 



2-N2 



(EQ.6) 



In Equation 6, pair (ij) represents a pixel of processed SAR data 204 within a block 
(which is a portion, A(ij) represents a corresponding in-phase or quadrature value of the pixel, 
and B(kiJc2) represents a corresponding DCT coefficient, where pair (kiMi) identifies the DCT 
5 coefficient in the DCT matrix. In the embodiment, a DCT coefficient is calculated over an eight 
by eight pixel block, i.e. Ni and N2 equal 8, alttiough other embodiments of the invention may 
use a different value for N. (The collection of DCT coefficients may be represented by an 8 by 8 
matrix.) 

Figure 7 shows a partitioning of complex unage data in order to obtain Discrete Cosine 
10 Transform (DCT) data in accordance with an embodiment of the invention. In the embodiment, 
SAR data 202 comprise the I component and the Q component, each component corresponding 
to a large (such as 1024 by 1024) data file 701. Transform module 203 partitions each file 701 
into a square (such as 8 by 8 block) (DCT matrix), e.g. blocks 703 and 705. Transform module 
203 processes each block (e.g. 703 and 705) in accordance with Equation 6. In order to process 
15 the entire data file 701, preprocessor 201 processes 128 partitions for both the I component and 
the Q component. 

Figure 8 shows apparatus 205 for quantizing Discrete Cosine Transform (DCT) data in 
accordance with an embodiment of the invention. Quantizer 205 comprises an adaptive table 
generation module 801, an adaptive table selection module 803, and a 

20 perform_data_quantization module 805. Adaptive table generation module 801 generates a new 
quantization conversion table (which contains a quantization matrix that is used for fiirther data 
compression as will be explained) for a new data type and stores the new quantization 
conversion table into a knowledge database 807 through an inter&ce 809 when functioning in a 
training mode but not during an operational mode. During the operational mode, transformed 

25 data 206 are processed by adaptive table selection module 803 and perfonn_data_quantization 
module 805. Depending upon the data type, as identified by adaptive control loop 251 from 
preprocessor 201, adaptive table selection module 803 selects an appropriate quantization 
conversion table, which comprises an 8 by 8 quantization matrix, firom knowledge database 807 
through an interface 811. If adaptive table selection module 803 cannot identify an appropriate 

30 quantization conversion table from adaptive control loop 251, module 803 selects a default 
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quantization conversion table. A quantization conversion table may correspond to different data 
types that are dependent upon factors including the type of radar, processing platform (which 
may affect the number of bits associated with SAR data 202), and topography that is associated 
withSAR data 202. 

5 Each element of a DCT matrix (e.g. matrix 703) is arithmetically divided by a 

corresponding element of the quantization matrix and rounded to an integer, thus providing 
quantized DCT data 208 (comprising a quantized I transform or a quantized Q transform). Each 
element of the quantization matrix is determined by statistics for the corresponding DCT 
coefficient in accordance with a specified maximum error (e.g. a root mean square error, a peak 

10 signal to noise ratio, and a byte by byte file comparison). (Figures 9 and 10 show statistics for 
the (1,1) and the (7,7) DCT coefficients, respectively.) The larger the value of an element of the 
quantization matrix, the greater flie corresponding step size (with less resolution). However, 
dividing an element of the DCT matrix by a larger number reduces the quantized value. If the 
quantized value is sufficiently reduced, the resulting value may be considered as being zero by 

1 5 encoder 207 if a specified maximum error (e.g. the root mean square error) is satisfied. 

In a variation of the embodiment, the quantization matrix may be determined by 
reducing a Measurement and Signature Intelligence (MASINT) product distortion. (In some 
cases, the reduction may correspond to a minimization of the distortion.) The distortion may be 
determined firom mterferometric SAR, coherent change detection (CCD), and polarimetry 

20 products. Interferometric SAR (IFSAR) is a comparison of two or more coherent SAR images 
collected at slightly different geometries. The process extracts phase differences caused by 
changes in elevation within the scene. IFSAR produces digital terrain elevation data suitable for 
use in providing terrain visualization products. (Products are generally referred as Digital 
Elevation Models (DEM).) These products are used in mapping and terrain visualization 

25 products. The advantage of IFSAR height determination is that is much more accurate than otheir 
methods, such as photo/radargrammetry methods that use only the intensity (magnitude) data, 
because phase is used and height determination is done with wavelength measurements which 
are very accurate (i.e. for commercial systems at C Band (5 GHz) approximately 5.3 cm)). 

Coherent Change Detection (CCD) is a technique involving the collection and 
30 comparison of a registered pair of coherent SAR images from approximately the same geometry 
collected at two different times (before and after an event). The phase information, not the 
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magnitude, is used to determine what has changed between the first and second collection. This 
can determine scene changes to the order of a wavelength (5.3 cm) and may denote ground 
changes/activity occurring between collections. 

Polarimetry products are generally collected using systems that can independently 
5 radiate and collect vertical and horizontal complex SAR data. This technique is accomphshed by 
alternately radiating vertical and horizontally polarized SAR pulses, receiving on both horizontal 
and vertical antennas, and saving the complex data fi:om each. The product formed is a unique 
target signature for objects with an associated complex polarized radar reflectance. This 
technique is used in many automatic target recognition systems (ATR). 

10 In a variation of the embodiment, each member of the quantization matrix (associated 

with a quantization conversion table) is determined by a heuristic process 1100 as shown in 
Figure 1 1. A quantization matrix for SAR data 202 may be determined by selecting an element 
of the quantization matrix in step 1103 and perturbing the value of the selected element in step 
1 105. In step 1 105, the selected element is incremented and decremented by incremental values. 

15 In step 1107, root means square errors (RSME) are calculated for different compression ratios. 
The selected value of the selected element is the value corresponding to a minimal root mean 
square error. If there are more elements in the quantization matrix to be processed, as 
determined in step 1109, the element indices (ij) are incremented in step 1111, and the next 
element is selected in step 1103. Steps 1105 and 1107 are repeated for the next element. The 

20 calculation of the quantization matrix is completed after all the elements are processed. 

In another variation of the embodiment, the quantization matrix is determined by the 
statistical characteristics of the DCT matrix, as was previously discussed. The quantization 
matrix is subsequently modified according to heuristic process 1 100. 

Transformed data 206 are quantized according to corresponding transform statistics that 
25 are associated with the DCT coefficients. DCT coefficients can be represented as departures 
firom a standard statistical distribution function (e.g., Laplacian, Gaussian, or Rayleigh). (A 
Laplacian function has a form of e""^', while a Gaussian function has a form of e''^.) Figure 9 
shows a representative histogram for a low order Discrete Cosine Transform (DCT) coefficient, 
DCT coefficient (1,1), in accordance with an embodiment of the invention. A number of 
30 observations 901 is shown in relation to corresponding bin- values 903. Actual data 905 is shown 
along with a Laplacian relation 907 and a Gaussian relation 909. Also, Figure 10 shows a 
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representative histogram for a high ordCT Discrete Cosine Transform (DCT) coefficient, DCT 
coefficient (7,7), in accordance with an emhodiment of the invention. A numher of observations 
1001 is shown in relation to corresponding bm values 1003. Actual data 1005 is shown along 
with a Laplacian relation 1007 and a Gaussian relation 1009. Analysis of the exemplary SAR 

5 data reveals a relationship with respect to the low order and high order DCT coefficients. By 
plotting the Lsqplacian and Gaussian functions and comparing the corresponding values with the 
DCT coefficient data of tiie exemplary SAR data, it is determined that low order tenns can be 
better represented by the Laplacian function, and the higher order temis can be better 
represented by the Gaussian function for typical SAR data. Qxxantization by quantizer 205 is 

10 designed by accounting for the complex SAR image DCT statistics as exemplified by Figures 9 
and 10. As the probability distribution becomes more focused about a zero value for a DCT 
coefficient, the less is the relative significance of the DCT coefficient witii respect to otiier DCT 
coefficients. Consequentiy, the corresponding entry in the quantization conversion table may be 
greater for the DCT coefficient. 

15 Oflier embodiments of the invention may utilize other transform types such as a Discrete 

Fourier Transform (DFT) or a discrete z-transform, both transforms being well known in the art. 
However, witii a selection of a different transform, the transform statistics may be different as 
reflected by the design of quantizer 205. . 

Quantized SAR data 208 are consequently processed by encoder 207 (e.g. a Hufftnan 
10 encoder). Each output 210 (comprising a compressed I component and a compressed Q 
component) of encoder 207 comprises an encoder pair (comprising a number of zeros that 
precede output 210 and a number of bits that represent a value of the corresponding DCT 
coefficient) and the value of the corresponding DCT coefficient (SAR data 208). Encoder 207 
may provide additioiud compression by removing a degree of redundancy that is associated with 
25 the encoder pair and SAR data 208 (in which firequentiy occurring data strings that are 
associated with the quantized DCT coefficients are replaced with shorter codes). Other 
embodiments of the invention may utilize other types of encoders such as Shannon Fano coding 
and Arithmetic coding. Encoder 207 provides encoded data 210 to post-processor 209. 

Post-processor 209 may further process encoded data 210 in order to format data 210 
30 into a format that is required for storing (that may be associated with archiving compressed data) 
or for transmitting compressed data 212 through a communications device. In the embodiment, 
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the communications device may be a radio frequency transmitter that transmits from a plane to a 
monitoring station, utilizing a radio data protocol as is known in the art In the embodiment, for 
example, post-processor 209 may format a data file (corresponding to a S AR image) into records 
that can be accormnodated by a storage device. Also, post-processor may include statistical 
5 information and the data type regarding SAR data 202, The statistical information and the data 
type may be used for decompressing compressed SAR data 212 at a subsequent time. 

Compressed data 212 may be subsequently decompressed by using apparatus that utilizes 
inverse operations corresponding to the operations that are provided by apparatus 200 in a 
reverse order. Figure 12 shows an apparatus 1200 for decompressing compressed data 1212 (that 

10 was compressed by apparatus 200 as shown in Figure 2) into a decompressed data 1202 in 
accordance with an embodiment of the invention, (Decompressed data 1202 approximates data 
202 within a specified maximum error.) An inverse post-processor 1209, a decoder 1107, an 
inverse quantizer 1205, an mverse transform module 1203, and an inverse preprocessor 1201 
correspond to post-processor 209, encoder 207, quantizer 205, transform module 203, and 

15 preprocessor 201, respectively. However, an inverse operation may not be able to exactly 
recover data because of quantization restraints. For example, quantizer 203 divides a DCT 
coefficient by a corresponding element in the quantization matrix (which is obtaioed from a 
quantization conversion table selected by module 803 from knowledge database 807) and 
rounded to an integer. The operation of rounding to an integer may cause information about the 

20 DCT coefficient to be lost. Consequently, the lost information cannot be recovered by mverse 
quantizer 1205 in determining the DCT coefficient 

In the embodiment, compressed data 212 may be compliant with National Imagery 
Transmission Format (NTTF) standards, in which header information about user-defined data 
(e.g. a quantization matrix) may be included. Thus, compressed data 212 may be compatible 
25 with processing software in accordance with Joint Photographic Experts Group (JPEG) 
compression standards. 

Other embodiments of the invention may compress and decompress data that are 
characterized by a different number of components (often referred as dimensions). Data that is 
characterized by more than one component (e.g. 2, 3, or more components) are often referred as 
30 multidimensional data. In such cases, preprocessor 201 may determine statistical characteristics 
associated with each of the components and map each of the components to bins in accordance 
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with the statistical characteristics. Transform module 203 transfonns each of the components 
according to a selected transform (e.g. a Fast Fourier Transform). Quantizer 205 subsequently 
quantizes each of the transformed components. 

Classical Electro-Optical (EO) based metrics, such as root mean square error (RMSE) 
5 and Peak Signal to Noise Ratio (PSNR), are useful for evaluating the magnitude imagery, but 
the EO-based metrics may not provide sufficient information about the phase data or the other 
derived products. EO-based metrics provide a necessary but not a sufficient condition for 
complex data compression fidelity. Useful magnitude imagery may also be available firom the 
compression process. The processes that generate phase data driven products such as 
10 interferometry, CCD and polarimetry may be included in the evaluations. Additional SAR data 
product metrics may be implemented to evaluate the phase information and any degradation of 
tiie products caused by compression. n 

An evaluation of compressed exemplary SAR data as processed by apparatus 200 
indicates that, with SAR data 202 being compressed at ratios greater than twenty to one, 
15 apparatus 200 may achieve near-lossless results for magnitude images and minimal degradation 
to phase information. 

As can be appreciated by one skilled in the art, a computer system with an associated 
computer-readable medium containing instructions for controUmg the computer system can be 
utilized to implement the exemplary embodiments that are disclosed herein. The computer 
20 system may include at least one computer such as a microprocessor, digital signal processor, and 
associated peripheral electronic circuitry. 
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LWe Claim: 

1. A method for compressing received data comprising a received In-phase (I) 
component and a received Quadrature (Q) component, tiie method comprising the steps of: 

(a) determining a first statistical characteristic associated with the received I 
component and a second statistical characteristic associated with the received Q 
component; 

(b) mapping the received I component to one of a first plurality of bins in 
accordance with the first statistical characteristic in order to form a converted I 
component, wherein a first resolution of the converted I component is reduced with 
respect to the received I component; and 

(c) mapping the received Q component to one of a second plurality of bins in 
accordance with the second statistical characteristic in order to form a converted Q 
component, wherein a second resolution of the converted Q component is reduced with 
respect to the received Q component. 

2. The method of Claim 1, wherein the first statistical characteristic and the second 
statistical characteristic are essentially the same. 

3. The method of Claim 2, wherein one of the first and second statistical 
characteristics is approximated by another of the first and second statistical characteristics. 

4. The method of Claim 1, wherein step (a) comprises the steps of: 

(i) detemiining the first statistical characteristic firom a first plurality of data 
points that are associated with the received I component; and 

(ii) detemaining the second statistical characteristic fi'om a second plurality of 
data points that are associated with the received Q component. 

5. The method of Claim 1 , further comprising the steps of: 

(d) transforming the received I component and the received Q component in 
order to modify the first and second statistical characteristics; and 
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(e) inverse transforming ttie converted I component and the converted 
received Q component. 

6. The method of Claim 5, wherein step (d) utilizes a Fast Fourier Transform (FFT) 
and step (e) utilizes an inverse Fast Fourier Transform (IFFT). 

7. The method of Claim 1, further comprising the steps of: 

(d) transforming the converted I component into a transformed I component 
and the converted Q component into a transformed Q component; and 

(e) quantizing the transformed I component into a quantized I transform and 
the transformed Q component into a quantized Q transform. 

8. The method of Claim 7, wherein step (d) utilizes a Discrete Cosine Transform 

(DCT). 

9. The method of Claim 7, further comprising the step of: 

(f) encoding the quantized I transform into a compressed I component and 
the quantized Q transform into a compressed Q component. 

10. The method of Claim 9, wherein the compressed I component and the 
compressed Q component are compliant to National Imagery Transmission Format (NITF) 
standards. 

11. The method of Claim 7, wherein step (e) utilizes a quantization conversion table 
that is determined by a statistical characterization of the transformed I and Q components. 

12. The method of Claim 7, wherein step (e) utilizes a quantization conversion table 
that is determined by a heuristic process. 

13. The method of Claim 7, wherein step (e) utilizes a quantization conversion table 
,that is determined by a statistical characterization of the transformed I and Q components and by 
a heuristic process. 
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14. The method of Claim 7, wherem step (e) utilizes a quantization conversion table 
that is detemiined by reducing a Measurement and Signature Intelligence (MASDSfT) product 
distortion. 

5 

15. The method of Claim 1, wherein the received data comprises synthetic aperture 
radar (SAR) data, 

16. The method of Claim 1, wherein bins of the first plurality of bins and the second 
10 plurality of bins are assigned in accordance with a distribution of data points. 

17. The method of Claim 1, wherein the first statistical characteristic comprises a 
first Probability Density Function (PDF) that is associated with the received I component and 
the second statistical characteristic comprises a second Probability Density Function (PDF) that 

15 is associated with the received Q component 

18. The method of Claim 17, wherein a fu:st width of a first bin is inversely related to 
a first value of the first probability density function, and wherein a second width of a second bin 
is inversely related to a second value of the second probability density fimction. 

20 

19. The method of Claim 1, wherem step (b) is based upon a first dynamic range that 
is associated with the received I component and step (c) is based upon a second dynamic range 
that is associated with the received Q component, 

25 20. A computer-readable medium having computer-executable instructions for 

performing the steps recited in Claim 1 . 

21. A computer-readable medium having computer-executable instructions for 
performing the steps recited in Claim 7. 

30 

22. A method for compressing received data comprising a received In-phase (I) 
component and a received (Quadrature (Q) component, the method comprismg the steps of: 
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(a) txansfonning the received I component and the received Q component 
into a magnitude component and a phase component; 

(b) selecting a first number of bits firom the magnitude component and a 
different second number of bits firom the phase component; and 

(c) forming a converted I component and a converted Q component fix)m the 
first number of bits and tiie different second number of bits, 

23. The method of Claim 22, further comprising the stq)S of: 

(d) transforming the converted I component into a transformed I component 
and the converted Q component into a transformed Q component; and 

(e) quantizing the transformed I component into a quantized I transform and 
the transformed Q component into a quantized Q transform. 

24. The method of Claim 23, further comprising the step of: 

(f) encoding the quantized I transform into a compressed I component and 
the quantized Q transform into a compressed Q component. 

25. A computer-readable medium having computer-executable instructions for 
performing the steps recited in Claim 22. 

26. A computer-readable medium having computer-executable instructions for 
performing the steps recited m Claim 23. 

27. A method for decompressing data in order to approximate original data, the 
original data comprising an original In-phase (I) component and an original Quadrature (Q) 
component, the method comprising the steps of: 

(a) obtaining a converted In-phase (I) component and a converted Quadrature 
(Q) component; 

(b) determining a first statistical characteristic that is associated with the 
original I component and a second statistical characteristic that is associated with the 
original Q component; 
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(c) convCTting the converted I component into a decompressed I component 
and the converted Q component into a deconq)ressed Q component in accordance with 
the first and second statistical characteristics, wherein the decompressed I component 
and the decompressed Q component comprise a greater number of bits than the 
converted I component and the converted Q component, respectively. 

28. The method of Claim 27, wherein step (c) comprises the steps of: 

(i) mapping the converted I component to the decompressed I component 
jfrom a first bin in accordance witii the first statistical characteristic; and 

(ii) mapping the converted Q component to the decompressed Q component 
fi^om a second bin in accordance with the second statistical characteristic. 

29. The method of Claim 28, wherein the first statistical characteristic comprises a 
first probability density fimction that is associated with the original I component, and wherein 
the second statistical characteristic comprises a second probability density fimction that is 
associated with the original Q component. 

30. The method of Claim 27, wherein the original data comprises synthetic aperture 
radar (SAR) data. 

3 1 . The method of Claim 27, wherein step (c) comprises the steps of: 

(i) converting the converted I component and the converted Q component 
into a magnitude component and a phase component; and 

(ii) forming the decompressed I component and the decompressed Q 
component from the magnitude component and tiie phase component 

32. The method of Claim 27, wherein step (a) comprises the steps of: 

(i) obtaining a quantized I transform and a quantized Q transform; 

(ii) inverse quantizing the quantized I transform into a transformed I 
component and the quantized Q transform into a transfonned Q component; and 

(iii) inverse transforming the transformed I component into the converted I 
component and the transformed Q component into the converted Q component. 
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33. The method of Claim 32, wherein step (i) comprises the steps of: 

(1) obtaining a compressed I component and a compressed Q component; and 

(2) decoding the compressed I component into the quantized I transform and 
the compressed Q component into the quantized Q transforai. 

34. A computer-readable medium having computer-executable instructions for 
performing the steps recited in Claim 27. 

35. A computer-readable medium having computer-executable instructions for 
performing the steps recited in Claim 32. 

36. An apparatus for compressing received data, the received data comprising a 
received In-phase Q) component and a received Quadrature (Q) component, the apparatus 
comprising: 

a preprocessor that deteraiines a first statistical characteristic associated with the 
received I component and a second statistical characteristic associated with the received 
Q component and that converts the received I component into a converted I component 
in accordance with the first statistical characteristic and the received Q component into a 
converted Q component in accordance with a second statistical characteristic, wherein a 
first resolution of the converted I component and a second resolution of the converted Q 
component are reduced with respect to the received I and Q components, respectively; 

a transform module that transforms the converted I component into a transformed 
I component and the converted Q component into a transformed Q component; and 

a quantizer that obtains an indicator firom the preprocessor about a data type of 
the received data, selects a quantization conversion table in accordance with the 
indicator, and quantizes the transformed I component into a quantized I transform and 
the transformed Q component into a quantized Q transform by utilizing the quantization 
conversion table. 

37. The apparatus of Claim 36, wherein the first statistical characteristic and the 
second statistical characteristic are essentially the same. 
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38. The apparatus of Claim 36, further comprising: 

an encoder that removes a degree of redundancy jfrom the quantized I transform 
to form a compressed I component and from the quantized Q transform to form a 
compressed Q transform, 

5 

39. The apparatus of Claim 36, wherein the first statistical characteristic comprises a 
first probability density function that is associated with the received I component and the second 
statistic characteristic comprises a second probability density fimction that is associated with the 
received Q component. 

10 

40. The apparatus of Claim 36, wherein the received data comprises synthetic 
aperture radar (SAR) data. 

41. An apparatus for decompressing data, comprising: 

15 an inverse quantizer that obtains a quantized I transform and a quantized Q 

transform, obtains an indicator about a data type that is associated with the data, selects 
a quantization conversion table in accordance with the indicator, and inverse quantizes 
the quantized I transform into a transformed I component and the quantized Q transform 
into a transformed Q component; 

20 an inverse transform module that inverse transforms the transformed I component 

into a converted I component and the transformed Q component into a converted Q 
component; and 

an inverse preprocessor that converts the converted I component into a 
decompressed I component and the converted Q component into a decompressed Q 
25 component, wherein the decompressed I component and the decompressed Q componait 

comprise a greater number of bits than the converted I component and the converted Q 
component, respectively. 

42. The apparatus of Claim 41, further comprising: 

30 a decoder that obtains a compressed I component and a compressed Q component 

and that converts the compressed I conq)onent into the quantized I transform and the 
compressed Q component into the quantized Q transform. 
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43. A method of compressing Synthetic Aperture Radar (SAR) data into compressed 
data, the SAR data comprising a received In-phase (T) component and a received Quadrature (Q) 
component, the method comprising the steps of: 

(a) determining a first statistical characteristic associated with the received I 
component and a second statistical characteristic associated with the received Q 
component; 

(b) mapping the received I component to a first bin in accordance with the 
first statistical characteristic in order to form a converted I component, wherem a first 
plurality of bins comprises the first bin, wherein a first resolution of the converted I 
component is reduced with respect to the received I component, and wherein the first 
plurality of bins is assigned according to a first distribution of received I component data 
points; 

(c) mapping the received Q component to a second bin in accordance with 
the second statistical characteristic in order to form a converted Q component, wherein a 
second plurality of bins comprises the second bin, wherein a first resolution of the 
converted Q component is reduced with respect to the received I component, and 
wherein the second plurality of bins is assigned according to a second distribution of 
received Q component data points; 

(d) transforming the converted I component into a transformed I component 
and the converted Q component into a transformed Q component; 

(e) quantizing the transformed I component into a quantized I transform and 
the transformed Q component into a quantized Q transform; and 

(f) encoding the quantized I transform into a compressed I component and 
the quantized Q transform into a compressed Q component. 

44. A method of compressing Synthetic Aperture Radar (SAR) data into compressed 
data, tiie SAR data comprising a received Ih-phase (T) component and a received Quadrature (Q) 
component, the method comprising the steps of: 

(a) transforming the received I component and the received Q component 
into a magnitude component and a phase component; 

(b) selecting a first number of bits from the magnitude component and a 
different second number of bits from the phase component; 
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(c) fonriing a converted I compon^t and a converted Q component from the 
first number of bits and the different second number of bits, wherein a first resolution of 
the converted I component is reduced with respect to the received I component and 
wherein a second resolution of the converted Q component is reduced with respect to the 
received Q component; 

(d) transforming the converted I component into a transfonned I component 
and the converted Q component into a transformed Q component; 

(e) quantizing the transformed I component into a quantized I transform and 
the transformed Q component into a quantized Q transform; and 

(f) encoding the quantized I transform into a compressed I component and 
the quantized Q transform into a compressed Q component 
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